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Abstract 

A construction of the Coulomb-Breit Hamiltonian for a pair of fermions, 
considered as a quantum two-body system, immersed in an arbitrary back- 
ground gravitational field described by Einstein's General Relativity is pre- 
sented. Working with Fermi normal coordinates for a freely falling observer 
in a spacetime region where there are no background sources and ignoring 
the gravitational back-reaction of the system, the effective Coulomb-Breit 
Hamiltonian is obtained starting from the S-matrix element corresponding 
to the one-photon exchange between the charged fermionic currents. The 
contributions due to retardation are considered up to order (f/c)^ and they 
are subsequently written as effective operators in the relativistic quantum 
mechanical Hilbert space of the system. The final Hamiltonian includes ef- 
fects linear in the curvature and up to order (v/c)'^. 

Key words: Effective potential, two-body atom, arbitrary gravitational 
background. 
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1. Introduction 

The connection between Einstein's General Relativity (GR), the dynam- 
ical theory of space-time, and Quantum Mechanics (QM) has a long history 
in physics and more than ever is nowadays a subject of intense research, 
both from theoretical and experimental perspectives, mainly in connection 
with possible deviations from GR and/or QM. Until now, experiments have 
already confirmed that inertia and Newtonian gravity affect quantum par- 
ticles, mainly electrons and neutrons, in ways that are fully consistent with 
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GR down to distances of the order of 10~^ cm. Gravitational-inertial fields 
leave their mark on particle wave functions in a variety of ways; particularly, 
they induce quantum phases that have been measured in some of the most 
renowned experiments on this topic [H, [s], 0, S, [oj . Recently, the purpose 
of tests of gravity has been mainly focused in determining the validity of 
the equivalence principle at the level of atoms, molecules, quasi-particles and 
antimatter. This translates into a need for higher precision tests of GR, and 
atomic physics together with quantum optics offer some of the most accurate 
results and promising scenarios. 

Most of the experimental tests of GR within the quantum realm arise 
from predictions on fully covariant wave equations where inertia and gravity 
appear as external classical fields, providing very valuable information on 
how Einstein's views carry through into the quantum world. However, many 
of the initially proposed and tested effects considered mainly one-particle sys- 
tems. This approach leaves aside the possibility of using the internal quan- 
tum structure of atoms or molecules as additional parameters. Nevertheless, 
this situation is drastically changing recently due to the fast development of 
matter wave interferometry 0, 11, IJ. 

Broadly speaking, previous work related to the description of atoms in 
a background gravitational field can be divided in three main branches, 
though with much interrelation among them: (i) the study of gravitational 
modifications to atomic spectra of mainly hydrogen-like atoms, which were 
considered as a quantum reduced mass revolving around some force center 

A two-body description of 
(ii) Restricted proofs of the 



lol ini, H, H H, [H, M, 113, [l8|Jl9lj2^^ 



the atom can be found in Refs. 22, 23 



24 



equivalence principle for test bodies made of classical electromagnetically in- 
teracting particles, within the PPN framework or within the THefi formalism 

The formulation 
iii) 



for static spherically symmetric spacetimes [25|, |26|, |27|, |28 



29, 3d 



of a quantum equivalence principle has been studied in Refs. 
Tests and studies of the gravitational red shift, where atomic clocks were de- 
scribed as quantum mechanical systems within the THefi or PPN formalisms 



31, 32, 33, 34, 35, 36, 37 



In this letter we consider the effect of a classical background gravita- 
tional field, described by GR, upon matter at the atomic level, where the 
use of QM becomes mandatory. We extend previous work in the following 
aspects: (i) We will consider the case of hydrogen-hke atoms described as 
two-body systems in the way first discussed in Ref. 23j. The inclusion of 



the position of the atom as a new degree of freedom may provide additional 
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couplings to probe the effects of gravity upon quantum objects in the region 
where tidal forces become important, (ii) The background gravitational field 



has no particular spatial symmetry, (iii) Following the steps of Ref. [38 
a quantum mechanical relativistic two-body effective potential in the back- 
ground gravitational field is produced, which takes into account retardation 
effects arising from the one-photon exchange between the currents. This 
constitutes the gravitational generalization of the relativistic Coulomb-Breit 
potential js^. 

Most of the general assumptions underlying this work are similar to those 
of Refs. , [isl, [lil, I23I, I25I, [26|: (i) The external gravitational field is well de- 



scribed by GR, so that it satisfies Einstein's field equations. Additionally, we 
ignore back-reaction effects, (ii) All the interacting particles are influenced 
by the background gravitational field at the quantum level, (iii) There exists 
an ideal freely falling observer in a region outside the sources of the back- 
ground gravitational field. This observer determines that, in his coordinate 
patch, the metric of the spacetime has the form 

9aP = VafS + QabfjcX^x" = Tj^p + , (l) 

where 

2 - Q 1 - Q 

QoaOb = —RoaOb {x ) , Qoabc = —-Roabc {x ) , Qabcd = —-R-abcd {x ) . (2) 

Here x'^ are three spatial coordinates, x° is the proper time of the observer 
and Raises is the projection of the background Riemann tensor on the or- 
thonormal tetrad carried by him, which in general depends on the proper 
time due to the motion of the observer. The weakness of the gravitational 
interaction induces very small corrections upon the observables, so that in 
our final results we preserve only quantities up to first order in Rap^s- The 
metric ([1]) is that of a freely falling observer using Fermi normal coordinates, 
which are appropriate to describe local experiments performed by inertial 



observers 40|]. (iv) The spatial extension and the time duration of events 



related to observations in the quantum system are very small compared with 



^Our metric has signature -2, 77"'^ = diag{—l, +1, +1, +1), e^^ is the tetrad such that 
ea^if^ ef = gf"", a = {0,1,2,3} = {0, a} are tetrad indices, while ^l = {0,1,2,3} = 
{0,a} are coordinate indices. Also eg'' e"i, = 5^. The symbol = denotes an equality 
including at most terms linear in the curvature. 
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the characteristic lengths and times of appreciable changes in the observer's 



Riemann tensor [15|. This allows well defined energy levels and the use of 
time-independent perturbation theory. Within this adiabatic approximation 
it is possible to ignore the x° dependence of the metric and of all the objects 
constructed from it: the time coordinate becomes just a fixed but arbitrary 



parameter 26|. (v) Conditions (i) and (iii) imply that the observer deter- 
mines that his Ricci tensor R^p is equal to zero. We also assume that during 
the measurements performed by the observer there are no particle creation 
effects due to the gravitational field. 

The coupled equations that describe the electromagnetic interaction be- 
tween spin 1/2 fermions in a background gravitational field, outside the grav- 
itational sources, are 

ich-^^D^i) - mc^-^ = q-f^-^A^, D^D^A" = qij^ip = J", (3) 

where we have chosen the Lorentz gauge D^A'^ = 0. Here = — 
^uj^^^Jbd is the covariant derivative written in general terms, cu"^^ = 
e"^ (D^e^'^) is the spin connection, and Jbd are the Lorentz group gener- 
ators in the corresponding representation for each field. 

2. The photon propagator 

The fundamental quantity required to construct an effective relativistic 
Hamiltonian describing the electromagnetic interaction between two charges 



is the curved space Feynman Green function G^iy'{x',x) for the photon |41 
We consider its Hadamard representation because it is valid in all space- 
times jisj . The differential equation for the electromagnetic Feynman Green 
function is 

□ C/^ix, X') = -K^fs (X, X') ^4 (X, X') , (4) 

where □ = g^^V is the parallel propagator bitensor between x' and 
x, and (^4 (x, x') is the invariant Dirac delta distribution for our observer. The 
solution to dl]), that vanishes at infinity for arbitrary time, is 
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The expressions for the parallel propagator (x, x') = e^y" (x) (x') are 

E 

f^oo = ^00 + I [^00 (x) + hoo {x')] , Koa = l [3/ioa (x) + /iQa (a^')] > (6) 

Here a{x,x') is the Synge world function, given by 

a (x, x') = ^r^^, (x - xT (x - x'Y + ^Q.au.A''' {x - x' f {x - x' f . (8) 



with 



A" 



- (x"^x^ + x'^x"') + - (x - x'y {x - x'y 

2 3 



(9) 



Note that expression contains terms of higher order than first in the 
curvature, arising from the denominator a + ie. We will keep the complete 
singular structure of the Green function because, in the following calculations, 
we will require expressions for the corresponding poles up to first order in 
hfj_i,. Let us remark that a{x, x') = o"(x, x', x° — x'°). 

3. The one-photon interaction 

In order to incorporate gravitational corrections into the electromagnetic 
interaction between two charged particles, we will use the S-Matrix method 



41[ as described in Ref.[38j, but generalized to a slightly curved space |42 
We need to evaluate the S-matrix element corresponding to the exchange of 
one photon between the two fermionic currents located at X2 and Xi, shown 
in Fig. (1), which is given by 

Sf^^ = ^ / (^2) G,, (X2, xi) Jf^ (xi) dV^ rfVi, (10) 

where G^^ {x2,Xi) is the electromagnetic Feynman Green function ([5]) and 
d^V = \/—g (x) d'^x is the invariant volume element at point x. The tran- 
sition current for each charge g„ has the form given in Eq. (|3]). As a con- 
sequence of the adiabatic approximation, the Dirac fermions which describe 



•^The orthonormal tetrad corresponding to the metric ([T]) is: ego = — 1 — 5^001 



^l^Qa, GaO — jhoa, Cab — Va 
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the charged particles can be represented by spinors with a definite positive 
frequency, in the coordinate patch of our observer. We can write the spinor 
field as 

^ (y) = e-^'^V^ (y) ^ith A;° > 0, (11) 

which satisfies the generalized Dirac equation given in The explicit form 
for each current is □ 



e'^n y , (12) 



where we have defined A^'^'' = fc^^"^ — k^f^\ Let us construct the two-particle 
states $jj (1,2) as the usual direct product of the spinorial wave function 
\l/^'^''(x) of each particle. We can then rewrite as a matrix element of an 
operator W^"'''^-' between the initial and final two-particle states : <l>j (2, 1) = 
^ff ^(x2) ® ^f^(xi) and $/ (2, 1) = ^(x2) ® ^^^^(xi), respectively, in such 
a way that 

% = / i^ff* (7f2)7?2)) (7fi)7?i)) d'V, d'V, ^ z6{n) ($/| >Ve|$.), 

(13) 

where n = Af^ + A^^ and 5{n) leads to energy conservation. Substituting 
(fT2|) in (fTOj) and with help of (fT3|) we identify the operator 



..,(2,1) ^ TkZik r (o^^2s}i!h. (^^^'^\ (iA\ 



where 



6;,,(X2, xi) = / / — --rfx^rfx? ^ r/,p / (2, 1) . (15) 

From expression f fl4|) for W'-^'"'"-* we will obtain the effective interaction op- 
erator, after evaluating the integrals in f lTSj) and after choosing an ordering 
prescription to complete the construction. To first order in the curvature, 
the poles of the Synge function are slightly displaced along the real line with 



■^Our conventions for the gamma matrices are: 7" denote fiat-space matrices, such, that 
{^a^^^l ^ 277"'5, (7°)^* = -7O, (75)^* ^ ^a. We have 7^(x) = ea^(a;)7". 
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respect to those of the flat expression. Because of this, we can take in f llSp 
the same integration contour used in the Minkoswki spacetime calculation, 
but inserting the modifled poles. The result of the contour integration is 



1(2,1) = -i6 (Q) 



bo 

tWj-: r 

e (i-<»oo) 



IpI (1 - ooo) 



[e (w) e™l^l + 9 i-w) e"™l''l] 



where 



IpI 



r 



1 + 



(16) 



(17) 



We have introduced the additional notation 



r = 5ab {,X2 - xiY {x2 - xi) , aoo = QoaObA' 



ab 



(18) 

bo = QoacbA''' - XiY , a = QcadbA'^' (^2 " X,)' (Xs - X^)" . (19) 

Here A"'^ is obtained from Eq. (|9]) by making the replacements x' ^ Xi, x ^ 
X2- The study of the symmetry properties under the exchange 1 -H- 2 in Eqs. 
(|T3l) and IHM leads to the property I (2, 1) = / (1, 2), which can be directly 
verifled in expression (JlGll . Next we expand both exponentials up to order 
w"^ obtaining 



-i6 (il) rjpp <j - 



2 ^"00 + ;:^ 



+ i 



wbo 
2r 



(20) 



The above expression still contains terms of order higher than flrst in the 
curvature because w = Wo + Wh- The flnal reduction to the desired order will 
be obtained after we rewrite the powers of w as functions of the one-particle 
Hamiltonian in the presence of a gravitational background 



(x) - ^h^opa + imc^ 



-Pa + aPc + ^hooPc 



1 



l--/ioo 7" + 7^c^/ioa7',(21) 



which is hermitian in the curved inner product 



(22) 
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Let us recall that in Eq. ( 1?T]) we have expanded all curved quantities up to 
first order in the curvature, which explains the contraction between tetrad 
and coordinate indexes that appears. This feature will frequently repeat 
itself in what follows. 



4. The effective relativistic interaction 

The effective interaction operator VVg is constructed in such a way that the 
w-dependent contributions in fl2U]) arise as a consequence of taking the matrix 
elements of We indicated in the right hand side of Eq. (fT3l) . The operator 
We will be a matrix operator acting on the quantum-mechanical Hilbert 
space of the system. To perform such construction we use the identification 
pa = —ihda in fl2ip . together with the eigenvalue equation for each spinor 
and Eq. (120]) . We start by substituting (1201) in (IT^ in such a way that, in 
the inner product (l22l) . we rewrite the transition element as 



Sji = i5{n){^f\d 



iw\ 1 f iw 



r (1 + F2 



l-^-.), (23) 



with the additional definitions 



O 




5c 



3^00 



1 



6 



C^OO o 



,(24) 
(25) 



Most applications will require to consider the low velocity limit of both par- 
ticles, which arises from the gamma matrices contributions to the operator O 
after taking the non-relativistic approximation. In order to keep track of the 
desired order in the approximation, it is convenient to expand O in powers 
of 

(26) 



A k 

A(„) 



which will produce terms of order f ^ 



-lln)l(n)^ 

j(^n)l^ corresponding limit. We 

consider this expansion up to second order. To preserve hermiticity in the 
inner product ( |22l) our construction must keep the expansion of all quantities 
up to first order in the curvature. In this way, the required expression for 
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O = Ea=o is given by 

1 



Oo-= 

6.1 = - 



I4x4 
1 



h',> + h 



0.= 



4J + C)(Af2) + Afi))-z6^-^, 



,(1) 



^{2) {hoc'^d)(i) + A('i) (/;,ocSj)(2) 



where the subindex A = 0,-1,-2 denotes the power to be achieved in 
the non-relativistic approximation]^ Substituting (127|) in (123|) . rewriting w 
as w {hc)~^, where u is associated to the difference of energy eigenvalues, we 
find 



with 



Wq = 6, 



e t^O (PC, 



6_i0c + 6o (^]f. 



0_ 



2<PC 



2^7 



'4x4 



Wc 



\2hcJ 



r(l + F2 



(28) 

(29) 
(30) 



We must now choose a prescription in order to convert each term of fl25]) into 
a hermitian operator acting on the quantum-mechanical Hilbert space of the 
system, in particular upon the wave functions in ( 128|) . Following Greiner 



38 , we take 



u^U = (2i) 



(31) 



where R and U are arbitrary operators, hermitian in the curved inner prod- 
uct, and "H*-") is the free Dirac Hamiltonian ( 12T]) for the n-th particle. The 
construction ( ETj) relies on the fact that both |<l>j), are eigenstates of 
the Hamiltonian for each particle. Before calculating the corresponding 
commutators we introduce some compact notation 



(32) 



(27) 



We have introduced the notation: h 



in) 
fiiy 



hi_,y (x„) , n = 1,2 and 



A a \b 



(n)c 
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X2 - Xi, 



z'' = x'i + x" 
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r 

n = -, 

r 

F2 = F2 



^00 



'00 ) 



1 



"-oo ^ "' 



(33) 
(34) 



Following (13T1) . we construct the operators 7^ : , Wf^ and arising 
from R : Wq, W-i and 1^-2, respectively. The final effective two-body inter- 
action due to one-photon exchange, which includes retardation effects up to 
second order, is 

^ Wf + >Vf ^ + >Vf . (35) 

By construction VVg is hermitian in the scalar product (122!) and we identify 
it with the interaction energy of the two charges, as seen by the inertia! 
observer. After a long and tedious calculation, keeping only linear terms in 
the curvature, we find 



-Of. 



1 

2r 



1 

' 2 Qoabd 



A 



ab 



cd 



1 + 



' 12 



\a^b xcd 



"-4x4 



Oac6 



AC A_ -y"- \ b 



( [2n"n^ - 35"^] - 4? 5"^) + ^Afa) [n"n'' - 35"^] (h x zS; 



n 



(1) 



1 I ~ '^'^ 

1 + floo - ^ 

,(1) 



i^2 Mfcc + 1 + nf,n, 



+ ((5^'^ + 2n'n'^) /i^^^ + 2n^e« + S;,, A^^. + (1^2 



(36) 



Equation (136!) constitutes the main result of this work. In the case of zero 
gravitational interaction, (136|) reduces to the flat space Coulomb plus Breit 
interaction. Contrary to the flat spacetime case. We contains terms of order 
arising from the contributions proportional to A^^-j. Those terms are 
essential in the hermiticity of the complete interacting term. The corrected 
effective potential exhibits two special features, which were absent in previ- 
ous discussions of the problem: (i) the coupling of spatial and spin degrees 
of freedom and (ii) the coupling of the standard center of mass and relative 



^The notation is: (h x S)" = e^'^^/iohSg. Since the results are to first order in the 
curvature, indices of quantities which are factors of h^^i, are lowered or raised by the flat- 
space metric. 
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coordinates. These effects, which also appear in the case of two charged par- 



ticles in the presence of an external inhomogeneous magnetic field [4^, may 
offer additional possibilities for gravitational testable effects at the quantum 
level. Using fl36|) we construct the relativistic Hamiltonian for two charges in 
a background gravitational field, including gravitational corrections due to 
one-photon exchange, as 

f|(l,2)=H(i)+H(2)+We, (37) 

where Eq. fl2ip gives each free Dirac Hamiltonian l-L. 

5. Final Comments 

Summarizing, for an observer attached to a freely falling frame outside the 
gravitational sources, we have constructed a quantum-mechanical two-body 
relativistic Hamiltonian describing the electromagnetic interaction between 
two charges in the presence of an arbitrary weak background gravitational 
field. Such Hamiltonian takes into account retardation effects up to order 
{v/cf' and it is hermitian in the curved scalar product (!22|) . The gravita- 
tional corrections are calculated only up to terms linear in the curvature. In 
order to extract physical consequences of the above Hamiltonian it is usu- 
ally necessary to take the non-relativistic approximation. In our case this 



will require to perform a two-body Foldy-Wouthuyzen transformation |45 
Even though this transformation can be performed in the full Hamiltonian 
fl37|) . the result is rather involved and does not shed much light unless one 



considers a particular situation of physical interest [46|. When dealing with 
applications, our general strategy will be to start by selecting and estimating 
the order of magnitude of the many competing terms that appear in (l37|l . 
This will allow us to isolate the dominant contributions of the gravitation- 
ally induced terms, which will be subsequently considered as perturbations 
of the resulting zeroth-order Hamiltonian. In many situations our observer 
would be interested in a gas of atoms at temperature T located at some 
distance Ir from the origin of the freely-falling frame. Five lengths are rele- 
vant and must be estimated: the local curvature radius D, the typical size 
of the coordinate patch Iq, the estimated position of the atom Ir, the inter- 
particle distance among the charges ao, and the de Broglie wavelength of 
the atom A^. The dimensionless parameters that measure the gravitational 
effects are 7' = {lo/Df , 7' = ilR/Df , 7o' = («o/^)' , 7' = {K/Df , 
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where we take Oq to be the Bohr radius of the atom. We certainly have 
7^,7q,7^ <<< 7^ << 1. Nevertheless, in comparison with the one-particle 
formulation of the problem, the presence of the additional parameter 7^, 
which is related to the location of the atom, allows us to probe regions where 
the tidal forces become relevant. Even considering its expected small value, 
this parameter could provide some amplifying factor leading to situations 
where observability might be enhanced. Given that we require to perform 
a non-relativistic approximation in our two-body relativistic Hamiltonian 
( 137|) . we must also consider additional dimensionless parameters related to 
the velocity of the atom rj = Vr ■ and to the velocity of the electron 
?7o ~ tte ~ 7.3 X 10^^. Under solar system conditions, for T ~ 10^ — 10^ K, 
we find rj ~ 10~^ — 10~^. Using these values of the parameters we can make 
further estimations: for astronomical observations within the solar system 
j2 y ^^^^ strong gravitational fields near the origin 7q might become the 
dominant parameter, and for very low temperature regimes 7^ could become 
relevant. Specific applications are beyond the scope of this letter and will be 



considered in forthcoming publications [46|, |47 
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Figure 1: The one-photon exchange diagram leading to the effective interaction potential. 
Fermions and photons include the effects of the gravitational background. 
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